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Decaying superheavy particles can be produced by Topological Defects or, in case they are quasi- 
stable, as relics from the early Universe. The decays of these particles can be the sources of observed 
Ultra High Energy Cosmic Rays (E ~ 10 10 - 10 12 GeV). The Topological Defects as the UHE CR 
sources are critically reviewed and cosmic necklaces and monopole-antiminopole pairs are identified 
as most plausible sources. The relic superheavy particles are shown to be clustering in the halo and 
their decays produce UHE CR without GZK cutoff. The Lightest Supersymmetric Particles with 
Ultra High Energies are naturally produced in the cascades accompanying the decays of superheavy 
qq ] particles. These particles are discussed as UHE carriers in the Universe. 
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i— I . I. INTRODUCTION 

a : 

The observation of cosmic ray particles with energies higher than 10 11 GeV Q gives a serious challenge to the 
known mechanisms of acceleration. The shock acceleration in different astrophysical objects typically gives maximal 
energy of accelerated protons less than (1 — 3) • 10 10 GeV [0. The unipolar induction can provide the maximal energy 
| 1 ■ 10 11 GeV only for the extreme values of the parameters^. Much attention has recently been given to acceleration 
by ultrarelativistic shocks jf|, ||. The particles here can gain a tremendous increase in energy, equal to T 2 , at a 
\ single reflection, where T is the Lorentz factor of the shock. However, it is known (see e.g. the simulation for pulsar 
relativistic wind in |^|) that particles entering the shock region are captured there or at least have a small probability 
to escape. 

Topological defects, TD, (for a review see |7|) can naturally produce particles of ultrahigh energies (UHE). The 
pioneering observation of this possibility was made by Hill, Schramm and Walker B (for a general analysis of TD as 
. UHE CR sources see || and for a review |l0|]). 

In many cases TD become unstable and decompose to constituent fields, superheavy gauge and Higgs bosons (X- 
particles), which then decay producing UHE CR. It could happen, for example, when two segments of ordinary string, 
or monopole and antimonopole touch each other, when electrical current in superconducting string reaches the critical 
^~ . value and in some other cases. 

In most cases the problem with UHE CR from TD is not the maximal energy, but the fluxes. One very general 
reason for the low fluxes consists in the large distance between TD. A dimension scale for this distance is the Hubble 
distance Hq . However, in some rather exceptional cases this dimensional scale is multiplied to a small dimensionless 
value r. If a distance between TD is larger than UHE proton attenuation length (due to the GZK effect jyj), then 
y\ ' the flux at UHE is typically exponential suppressed. 

Ordinary cosmic strings can produce particles when a loop annihilate into double line |l2| . The produced UHE CR 
flux is strongly reduced due to the fact that a loop oscillates, and in the process of a collapse the two incoming parts 
of a loop touch each other in one point producing thus the smaller loops, instead of two-line annihilation. However, 
this idea was recently revived due to recent work It is argued there that the energy loss of the long strings is 
dominated by production of very small loops with the size down to the width of a string, which immediately annihilate 
into superheavy particles. A problem with this scenario is too large distance between strings (of order of the Hubble 
distance). For a distance between an observer and a string being the same, the observed spectrum of UHE CR has 
an exponential cutoff at energy E ~ 3 ■ I0 19 eV . 

Superheavy particles can be also produced when two segments of string come into close contact, as in cusp events 
fwj. This process was studied later by Gill and Kibble [O], and they concluded that the resulting cosmic ray flux 
is far too small. An interesting possibility suggested by Brandenberger [ H is the cusp "evaporation" on cosmic 
strings. When the distance between two segments of the cusp becomes of the order of the string width, the cusp 
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may "annihilate" turning into high energy particles, which are boosted by a very large Lorentz factor of the cusp p4| . 
However, the resulting UHE CR flux is considerably smaller than one observed fl6| |. 

Superconducting strings flTj] appear to be much better suited for particle production. Moving through cosmic 
magnetic fields, such strings develop electric currents and copiously produce charged heavy particles when the current 
reaches certain critical value. The CR flux produced by superconducting strings is affected by some model-dependent 
string parameters and by the history and spatial distribution of cosmic magnetic fields. Models considered so far 
failed to account for the observed flux J[8| . 

Monopole-antimonopole pairs (MM) can form bound states and eventually annihilate into UHE particles [ ^9[ , 
. For an appropriate choice of the monopole density um, this model is consistent with observations; however, the 
required (low) value of hm implies fine-tuning. In the first phase transition G — > H x U(l) in the early Universe the 
monopoles are produced with too high density. It must then be diluted by inflation to very low density, precisely 
tuned to the observed UHE CR flux. 

Monopole- string network can be formed in the early Universe in the sequence of symmetry breaking 

G -► H x U(l) -> H x Z N . (I) 

For N > 3 an infinite network of monopoles connected by strings is formed. The magnetic fluxes of monopoles in the 
network are channeled into into the strings that connect them. The monopoles typically have additional unconfined 
magnetic and chromo-magnetic charges. When strings shrink the monopoles are pulled by them and are accelerated. 
The accelerated monopoles produce extremely high energy gluons, which then fragment into UHE hadrons pl| . The 
produced flux is too small to explain UHE CR observation . 

Cosmic necklaces are TD which are formed in a sequence of symmetry breaking given by Eq.(|) when N = 2. 
The first phase transition produces monopoles, and at the second phase transition each monopole gets attached to 
two strings, with its magnetic flux channeled along the strings. The resulting necklaces resemble "ordinary" cosmic 
strings with monopoles playing the role of beads. Necklaces can evolve in such way that a distance between monopoles 
diminishes and in the end all monopoles annihilate with the neighboring antimonopoles (2^] . The produced UHE CR 
flux [^3| is close to the observed one and the shape of the spectrum resembles one observed. The distance between 
necklaces can be much smaller than attenuation length of UHE protons. 

Superheavy relic particles can be sources of UHE CR |24|,|25| . In this scenario Cold Dark Matter (CDM) have a small 
admixture of long-lived superheavy particles. These particles must be heavy, rax > 10 12 GeV, long-lived tx > to, 
where to is the age of the Universe, and weakly interacting. The required life-time can be provided if this particle 
has (almost) conserved quantum number broken very weakly due to warmhole |2f| or instanton effects. Several 
mechanisms for production of such particles in the early Universe were identified. Like other forms of non-dissipative 
CDM , X-particles must accumulate in the halo of our Galaxy |25| and thus they produce UHE CR without GZK 
cutoff and without appreciable anisotropy. 

The UHE carriers produced at the decay of superheavy relic particles or from TD, can be be nucleons, photons 
and neutrinos or neutralinos |2§| l . Production of neutralinos occurs in particle cascade, which originates at the decay 
of superheavy X-particle in close analogy to QCD cascade ]2q| . Though flux of UHE neutralino is of the same order 
as neutrino flux, its detection is more problematic because of smaller cross-section. The other particles discussed as 
carrier of UHE signal are gluino ]2^-]28[|, in case it is Lightest Supersymmetric Particle (LSP), and heavy monopole 

us 

In this paper I will present the results obtained in our joint works with M.Kachelriess and A. Vilenkin about necklaces 
and superheavy relic particles as possible sources of UHE CR. Neutralino and gluino as the carriers of UHE signal 
will be also shortly discussed. 

II. NECKLACES 

Necklaces produced in a sequence of symmetry breaking G — » H X £7(1) H x Z2 form the infinite necklaces 
having the shape or random walks and a distribution of closed loops. Each monopole in a necklace is attached to two 
strings. 

The monopole mass m and the string tension are determined by the corresponding symmetry breaking scales, rj s 
and r] m (rj m > n s ): m ~ 4nri m /e, fi ~ ^vt- Here, e is the gauge coupling. The mass per unit length of string is 
equal to its tension, //. Each string attached to a monopole pulls it with a force F ~ /1 in the direction of the string. 
The monopole radius S m and the string thickness S s are typically of the order S m ~ (&q m ) , S s ~ (er] s ) . 

An important quantity for the necklace evolution is the dimensionless ratio 
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r = m/ fid, 



(2) 



We expect the necklaces to evolve in a scaling regime. If £ is the characteristic length scale of the network, equal 
to the typical separation of long strings and to their characteristic curvature radius, then the force per unit length of 
string is / ~ and the acceleration is a ~ (r + 1) We assume that £ changes on a Hubble time scale ~ t. 

Then the typical distance travelled by long strings in time t should be ~ £, so that the strings have enough time to 
intercommute in a Hubble time. This gives at 2 ~ £, or 

Z~( r + 1)-VH. (3) 

The typical string velocity is v ~ (r + l) -1 / 2 . 

It is argued in Ref.( j23|) that r(t) is driven towards large value r > 1. However, for r > 10 6 the characteristic 
velocity of the network falls down below the virial velocity, and the necklaces will be trapped by gravitational clustering 
of the matter. This may change dramatically the evolution of network. One possible interesting effect for UHE CR 
can be enhancement of necklace space density within Local Supercluster - a desirable effect as far absence of the GZK 
cutoff is concerned. However, we restrict our consideration by the case r < 10 6 . The distance between necklaces is 
still small enough, £ > 3 Mpc, to assume their uniform distribution, when calculating the UHE CR flux. 

Self-intersections of long necklaces result in copious production of closed loops. For r > 1 the motion of loops is 
not periodic, so loop self-intersections should be frequent and their fragmentation into smaller loops very efficient. A 
loop of size I typically disintegrates on a timescale r ~ r~ x / 2 l. All monopoles trapped in the loop must, of course, 
annihilate in the end. 

Annihilating MM pairs decay into Higgs and gauge bosons, which we shall refer to collectively as X-particles. The 
rate of X-particle production is easy to estimate if we note that infinite necklaces lose a substantial fraction of their 
length to closed loops in a Hubble time. The string length per unit volume is ~ £~ 2 , and the monopolc rest energy 
released per unit volume per unit time is rji/^ 2 t. Hence, we can write 

h x ~r 2 n/(t 3 m x ), (4) 

where nix ~ e?y m is the X-particle mass. 

X-particles emitted by annihilating monopoles decay into hadrons, photons and neutrinos, the latter two components 
are produced through decays of pions. 

The diffuse flux of ultra-high energy protons can be evaluated as 

ch x f ta ^ TI7 _ s dE g (E,t) 



*>W = 4^J dtW N (mx, Xg )^^ (5) 

where dnx/dt is given by Eq.(|J), E is an energy of proton at observation and E g (E, t) is its energy at generation at 
cosmological epoch t, x g — E g j E and Wximx, x) is the fragmentation function of X-particle into nucleons of energy 
E = xmx- The value of dE g ldE can be calculated from the energy losses of a proton on microwave background 
radiation (e.g. see ||). In Eq.(J5|) the recoil protons are taken into account, while in Ref.( |23jl ) their contribution was 
neglected. 

The fragmentation function Wximx, x) is calculated using the decay of X-particle into QCD partons (quark, gluons 
and their supersymmetric partners) with the consequent development of the parton cascade. The cascade in this case 
is identical to one initiated by e + e _ -annihilation. We have used the fragmentation function in the gaussian form as 
obtained in MLLA approximation in [Q and |35|| . 

In our calculations the UHE proton flux is fully determined by only two parameters, r 2 \x and mx- The former is 
restricted by low energy diffuse gamma-radiation. It results from e-m cascades initiated by high energy photons and 
electrons produced in the decays of X-particles. 

The cascade energy density predicted in our model is 

if r 2 u t d ± 1 = 3 -fr*» 

r M io * 3 (i + ^) 4 " ± t 2 



u cas = ^j^r ii \ — = -f n r -2, (6) 



where to is the age of the Universe (here and below we use h — 0.75), z is the redshift and f„ ~ 1 is the fraction of 
energy transferred to pions. In Eq.(0) we took into account that half of the energy of pions is transferred to photons 
and electrons. The observational bound on the cascade density, for the kind of sources we are considering here, is |36| ] 
u C as ^ 10~ 5 eV/cm 3 . This gives a bound on the parameter r 2 //. 
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In numerical calculations we used r 2 p = 1 x 10 28 GeV 2 , which results in ui cas = 5.6 • 10~ 6 eV/cm 3 , somewhat 
below the observational limit. Now we are left with one free parameter, mx, which we fix at 1 • 10 14 GeV. The 
maximum energy of protons is then E max ~ 10 13 GeV. The calculated proton flux is presented in Fig.l, together 
with a summary of observational data taken from ref . p7| . 

Let us now turn to the calculations of UHE gamma-ray flux from the decays of X-particles. The dominant channel 
is given by the decays of neutral pions. The flux can be readily calculated as 

I 7 (E) = ±-n x \ 7 (E)N 7 (E), (7) 

where fix is given by Eq.(Q), Xj(E) is the absorption length of a photon with energy E due to e + e~ pair production 
on background radiation and N 7 (E) is the number of photons with energy E produced per one decay of X-particle. 
The latter is given by 

2 f 1 clr 

N y (E) = / ^ Wn o(m x ,x) (8) 

m X J E /m x X 

where W 7r o(mx, x) is the fragmentation function of X-particles into 7r° pions. 

At energy E > 1 • 10 10 GeV the dominant contribution to the gamma-ray absorption comes from the radio back- 
ground. The significance of this process was first noticed in [Q(see also book S). New calculations for this absorption 
were recently done ||^] . We have used the absorption lengths from this work. 

When evaluating the flux| at E > MO 10 GeV we neglected cascading of a primary photon, because pair production 
and inverse compton scattering occur at these energies on radio background, and thus at each collision the energy 
of a cascade particle is halved. Moreover, assuming an intergalactic magnetic field H > 1 • 10~ 9 , the secondary 
electrons and positrons loose their energy mainly due to synchrotron radiation and the emitted photons escape from 
the considered energy interval [|Io| . 

The calculated flux of gamma radiation is presented in Fig.l by the curve labelled 7. One can see that at E ~ 
1 • 10 11 GeV the gamma ray flux is considerably lower than that of protons. This is mainly due to the difference in 
the attenuation lengths for protons (110 Mpc) and photons (2.6 Mpc [[39) and 2.2 Mpc fl38f ). At higher energy the 
attenuation length for protons dramatically decreases (13.4 Mpc at E = 1 • 10 12 GeV) and the fluxes of protons and 
photons become comparable. 

A requirement for the models explaining the observed UHE events is that the distance between sources must 
be smaller than the attenuation length. Otherwise the flux at the corresponding energy would be exponentially 
suppressed. This imposes a severe constraint on the possible sources. For example, in the case of protons with energy 
E ~ (2 — 3) • 10 11 GeV the proton attenuation length is 19 Mpc . If protons propagate rectilinearly, there should be 
several sources inside this radius; otherwise all particles would arrive from the same direction. If particles are strongly 
deflected in extragalactic magnetic fields, the distance to the source should be even smaller. Therefore, the sources 
of the observed events at the highest energy must be at a distance R < 15 Mpc in the case or protons. 

In our model the distance between sources, given by Eq.(||), satisfies this condition for r > 3 • 10 4 . This is in contrast 
to other potential sources, including supeconducting cosmic strings and powerful astronomical sources such as AGN, 
for which this condition imposes severe restrictions. 

The difficulty is even more pronounced in the case of UHE photons. These particles propagate rectilinearly and 
their absorption length is shorter: 2 — 4 Mpc at E ~ 3 • 10 11 GeV. It is rather unrealistic to expect several powerful 
astronomical sources at such short distances. This condition is very restrictive for topological defects as well. The 
necklace model is rather exceptional regarding this aspect. 



III. UHE CR FROM RELIC QUASISTABLE PARTICLES 

This possibility was recognized recently in Refs( pjL^Sj ). 

Our main assumption is that Cold Dark Matter (CDM) has a small admixture of long-lived supermassive X- 
particles. Since, apart from very small scales, fluctuations grow identically in all components of CDM, the fraction of 
X-particles, £x, is expected to be the same in all structures. In particular, £x is the same in the halo of our Galaxy 
and in the extragalactic space. Thus the halo density of A-particles is enhanced in comparison with the extragalactic 
density. The decays of these particles produce UHE CR, whose flux is dominated by the halo component, and therefore 
has no GZK cutoff. Moreover, the potentially dangerous e-m cascade radiation is suppressed. 
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First, we address the elementary-particle and cosmological aspects of a superheavy long-living particle. Can the 
relic density of such particles be as high as required by observations of UHE CR? And can they have a lifetime 
comparable or larger than the age of the Universe? 

Let us assume that A-particle is a neutral fermion which belongs to a representation of the SU(2) x U(l) group. 
We assume also that the stability of A-particles is protected by a discrete symmetry which is the remnant of a gauge 
symmetry and is respected by all interactions except quantum gravity through wormholc effects. In other words, our 
particle is very similar to a very heavy neutralino with a conserved quantum number, R' , being the direct analogue 
of i?-parity (see Q and the references therein). Thus, one can assume that the decay of A-particle occurs due 
to dimension 5 operators, inversely proportional to the Planck mass mpi and additionally suppressed by a factor 
exp(— 5), where S is the action of a wormhole which absorbs i?'-charge. As an example one can consider a term 

1 -§z^0exp(-S), (9) 



mpi 



where ^ describes X-particle, and is a SU(2) scalar with vacuum expectation value vew = 250 GeV. After 
spontaneous symmetry breaking the term (S) results in the mixing of A-particle and neutrino, and the lifetime due 
to X^v + q + q, e.g., is given by 



{G F v z EW ) 2 m\ 



where Gf is the Fermi constant. The lifetime tx > to for A-particle with mx > 10 13 GeV needs S > 44. This value 
is within the range of the allowed values as discussed in Ref. [ pl| . 

Let us now turn to the cosmological production of A-particles with mx > 10 13 GeV. Several mechanisms were 
identified in p5|, including thermal production at the reheating stage, production through the decay of inflaton field 
at the end of the "pre- heating" period following inflation, and through the decay of hybrid topological defects, such 
as monopoles connected by strings or walls bounded by strings. 

For the thermal production, temperatures comparable to mx are needed. In the case of a heavy decaying gravitino, 
the reheating temperature Tr (which is the highest temperature relevant for our problem) is severely limited to value 
below 10 8 — 10 10 GeV, depending on the gravitino mass (see Ref. |}2| and references therein). On the other hand, 
in models with dynamically broken supersymmetry, the lightest supersymmetric particle is the gravitino. Gravitinos 
with mass m 3 / 2 < 1 keV interact relatively strongly with the thermal bath, thus decoupling relatively late, and can 
be the CDM particle |3j| . In this scenario all phenomenological constraints on Tr (including the decay of the second 
lightest supersymmetric particle) disappear and one can assume Tr ~ 10 11 — 10 12 GeV. In this range of temperatures, 
A-particles are not in thermal equilibrium. If Tr < mx , the density nx of A-particles produced during the reheating 
phase at time tR due to a + a ^ X + X is easily estimated as 

n x {tR) ~ N a n 2 a axt R ex.p(-2mx /T R ), (11) 

where N a is the number of flavors which participate in the production of X-particles, n a is the density of a-particles 
and ax is the production cross-section. The density of A-particles at the present epoch can be found by the standard 
procedure of calculating the ratio nx/s, where s is the entropy density. Then for mx — 1 ■ 10 13 GeV and £x m the 
wide range of values 10~ 8 — 10~ 4 , the required reheating temperature is Tr ~ 3 • 10 11 GeV. 

In the second scenario mentioned above, non-equilibrium inflaton decay, A-particles are usually overproduced and 
a second period of inflation is needed to suppress their density. 

Finally, A-particles could be produced by TD such as strings or textures. Particle production occurs at string 
intersections or in collapsing texture knots. The evolution of defects is scale invariant, and roughly a constant number 
of particles v is produced per horizon volume t 3 per Hubble time t. {v ~ 1 for textures and y> 1 for strings.) The 
main contribution to to the X-particle density is given by the earliest epoch, soon after defect formation, and we find 
£x ~ 10~ e v(mx /10 13 GeV)(Tf /10 10 GeV) 3 , where Tf is the defect formation temperature. Defects of energy scale 
V m x could be formed at a phase transition at or slightly before the end of inflation. In the former case, Tf ~ Tr , 
while in the latter case defects should be considered as " formed" when their typical separation becomes smaller than t 
(hence Tf < Tr). It should be noted that early evolution of defects may be affected by friction; our estimate of £x will 
then have to be modified. X particles can also be produced by hybrid topological defects: monopoles connected by 
strings or walls bound by strings. The required values of Tlx I s can be obtained for a wide range of defect parameters. 

The decays of A-particles result in the production of nucleons with a spectrum W]\[(mx, x), where mx is the mass 
of the X-particle and x = E/mx- The flux of nucleons (j),p,n,n) from the halo and extragalactic space can be 
calculated as 
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I i N (E) = ^-^R i —W N (m x ,x), (12) 

where index i runs through h (halo) and ex (extragalactic), Ri is the size of the halo Rh, or the attenuation length 
of UHE protons due to their collisions with microwave photons, X P (E), for the halo case and extragalactic case, 
respectively. We shall assume mxn x — C^Pcdm an d m xn e £ = £x^CDM/?cr, where £x describes the fraction of 
X-particles in CDM, Qcdm is the CDM density in units of the critical density p cr , and Pcdm ~ 0-3 GeV/cv 3 
is the CDM density in the halo. We shall use the following values for these parameters: a large DM halo with 
Rh = 100 kpc (a smaller halo with Rh = 50 kpc is possible, too), SlcDAih 2 — 0.2, the mass of V-particle in the 
range 10 13 GeV < mx < 10 16 GeV, the fraction of V-particles £x <C 1 and tx 3> to, where to is the age of the 
Universe. The two last parameters are convolved in the flux calculations in a single parameter rx — £,xto/ T x- For 
W]\[(mx, x) we shall use like in the previous section the QCD fragmentation function in MLLA approximation. For 
the attenuation length of UHE protons due to their interactions with microwave photons, we use the values given in 
the book §. 

The high energy photon flux is produced mainly due to decays of neutral pions and can be calculated for the halo 
case as 



I h ,{E) = ^^R h N 1 {El (i:S) 



where N 1 (E) is the number of photons with energy E produced per decay of one V-particle, which is given by Eq.(|^) 
For the calculation of the extragalactic gamma-ray flux, it is enough to replace the size of the halo, Rh, by the 
absorption length of a photon, X 7 (E). The main photon absorption process is e + e~-production on background 
radiation and, at E > 1 • 10 10 GeV, on the radio background. The neutrino flux calculation is similar. 
Before discussing the obtained results, we consider the astrophysical constraints. 

The most stringent constraint comes from electromagnetic cascade radiation, discussed in the previous section. In 
the present case this constraint is weaker, because the low-energy extragalactic nucleon flux is ~ 4 times smaller than 
that one from the Galactic halo (see Fig. 2). Thus the cascade radiation is suppressed by the same factor. 

The cascade energy density calculated by integration over cosmological epochs (with the dominant contribution 
given by the present epoch z = 0) yields in our case 

1 ^CDMPcr 1n2 , 3 

Wcas = -rx — — = 6.3 • 10 rxf-K eV/cm . (14) 

5 -Ho^o 

To fit the UHE CR observational data by nucleons from halo, we need rx = 5 • 10~ n . Thus the cascade energy 
density is ui ca , s = 3.2 • 10 _8 / w eV/cm 3 , well below the observational bound. 

Let us now discuss the obtained results. The fluxes shown in Fig. 2 are obtained for Rh — 100 kpc, mx = 1T0 13 GeV 
and rx = ^xto/rx = 5 ■ 10~ n . This ratio rx allows very small £x and tx > to- The fluxes near the maximum energy 
E max — 5 ■ 10 12 GeV were only roughly estimated (dotted lines on the graph). 

It is easy to verify that the extragalactic nucleon flux at E < 3 • 10 9 GeV is suppressed by a factor ~ 4 and by 
a much larger factor at higher energies due to energy losses. The flux of extragalactic photons is suppressed even 
stronger, because the attenuation length for photons (due to absorption on radio-radiation) is much smaller than for 
nucleons (see Ref. pfij]). This flux is not shown in the graph. The flux of high energy gamma-radiation from the halo 
is by a factor 7 higher than that of nucleons and the neutrino flux, given in the Fig. 2 as the sum of the dominant halo 
component and subdominant extragalactic one, is twice higher than the gamma-ray flux. 

The spectrum of the observed EAS is formed due to fluxes of gamma-rays and nucleons. The gamma-ray contribution 
to this spectrum is rather complicated. In contrast to low energies, the photon-induced showers at E > 10 9 GeV 
have the low-energy muon component as abundant as that for nucleon- induced showers [0]. However, the shower 
production by the photons is suppressed by the LPM effect [H and by absorption in geomagnetic field (for recent 
calculations and discussion see |4l],|^| and references therein). 

We wish to note that the excess of the gamma-ray flux over the nucleon flux from the halo is an unavoidable feature 
of this model. It follows from the more effective production of pions than nucleons in the QCD cascades from the 
decay of V-particle. 

The signature of our model might be the signal from the Virgo cluster. The virial mass of the Virgo cluster is 
Afyirgo ~ 1 • 10 15 M Q and the distance to it R = 20 Mpc. If UHE protons (and antiprotons) propagate rectilinearly 
from this source (which could be the case for E p ~ 10 11 — 10 12 GeV), their flux is given by 

Vn g o = rx _M^ WN[m ^ x) _ (15) 

^' zqR m x 
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The ratio of this flux to the diffuse flux from the half hemisphere is 6.4- 10 3 . This signature becomes less pronounced 
at smaller energies, when protons can be strongly deflected by intergalactic magnetic fields. 

IV. LSP IS UHE CARRIER 

LSP is the Lightest Supersymmetric Particle. It can be stable if R-parity is strictly conserved or unstable if R-parity 
is violated. To be able to reach the Earth from most remote regions in the Universe, the LSP must have lifetime 
longer than tlsp ^ io/T, where to is the age of the Universe and T = E/tiilsp is the Lorentz-factor of the LSP. In 
case rriLSP ~ 100 GeV, tlsp > 1 yr. 

Theoretically the best motivated candidates for LSP are the neutralino and gravitino. We shall not consider the 
latter, because it is practically undetectable as UHE particle. 

In all elaborated SUSY models the gluino is not the LSP. Only, if the dimension-three SUSY breaking terms are set 
to zero by hand, gluino with mass rrig = 0(1 GeV) can be the LSP There is some controversy if the low-mass 



window 1 GeV < m~ g < 4 GeV for the gluino is still allowed 46 47|] . Nevertheless, we shall study the production of 
high-energy gluinos and their interaction with matter being inspired by the recent suggestion p(J (see also H?]]), that 
the atmospheric showers observed at the highest energies can be produced by colorless hadrons containing gluinos. 
We shall refer to any of such hadron as g-hadron. Light gluinos as UHE particles with energy E > 10 16 eV were 
considered in some detail in the literature in connection with Cyg X-3 M,E9|. Additionally, we consider heavy gluinos 



with m~ g > 150 GeV |7| 



UHE LSP are most naturally produced at the decays of unstable superheavy particles, either from TD or as the 
relic ones [ p8| . 

The QCD parton cascade is not a unique cascade process. A cascade multiplication of partons at the decay of 
superheavy particle appears whenever a probability of production of extra parton has the terms alnQ 2 or a In 2 Q 2 , 
where Q is a maximum of parton transverse momentum, i.e. mx in our case. Regardless of smallness of a, the cascade 
develops as long as a\nQ 2 > 1. Therefore, for extremely large Q 2 we are interested in, a cascade develops due to 
parton multiplication through SU(2) x U(l) interactions as well. Like in QCD, the account of diagrams with alnQ 2 
gives the Leading Logarithm Approximation to the cascade fragmentation function. 

For each next generation of cascade particles the virtuality of partons q 2 diminishes. When q 2 3> rn 2 SUSY > where 
msuSY is a typical mass of supersymmetric particles, the number of supersymmetric partons in the cascade is the 
same as their ordinary partners. At q 2 < m 2 SUSY the supersymmetric particles are not produced any more and the 
remaining particles decay producing the LSP. In Ref.( pl|]) a simple Monte Carlo simulation for SUSU cascading 
was performed and the spectrum of emitted LSP was calculated. LSP take away a considerable fraction of the total 
energy (~ 40%). 

The fluxes of UHE LSP are shown in Fig. 3 for the case of their production in cosmic necklaces (see section II) . When 
the LSP is neutralino, the flux is somewhat lower than neutrino flux. The neutralino-nucleon cross- section, <j x n, is also 
smaller than that for neutrino. For the theoretically favorable masses of supersymmetric particles, <j x n ~ 10 -34 cm 2 
at extremely high energies. If the the masses of squarks are near their experimental bound, ~ 180 GeV , the 

cross-section is 60 times higher. 

Gluino as the LSP is another phcnomcnological option. Let us discuss shortly the status of the gluino as LSP. 

In all elaborated SUSY models the gluino is not LSP, and this possibility is considered on purely phenomenological 
basis. Accelerator experiments give the lower limit on the gluino mass as rrig > 150 GeV ^6|. The upper limit of 
the gluino mass is given by cosmological and astrophysical constraints, as was recently discussed in [ p7f . In this work 
it was shown that if the gluino provides the dark matter observed in our galaxy, the signal from gluino annihilation 
and the abundance of anomalous heavy nuclei is too high. Since we are not interested in the case when gluino is DM 
particle, we can use these arguments to obtain an upper limit for the gluino mass. Calculating the relic density of 
gluinos (similar as in f2^| ) and using the condition fig <C ricDM, we obtained wig <C 9 TeV. 

Now we come to very interesting argument against existence of a light stable or quasistable gluino |5(J . It is plausible 
that the glueballino (gg) is the lightest hadronic state of gluino j4q^9]. However, gluebarino, i.e. the bound state of 
gluino and three quarks, is almost stable because baryon number is extremely weakly violated. In Ref. fl50| it is argued 
that the lightest gluebarino is the neutral state (guud). These charged gluebarinos are produced by cosmic rays in 
the earth atmosphere |50| , and light gluino as LSP is excluded by the search for heavy hydrogen or by proton decay 



experiments (in case of quasistable gluino) . In the case that the lightest gluebarino is neutral, see 45 1 , the arguments 
of ]50[ | still work if a neutral gluebarino forms a bound state with the nuclei. Thus, a light gluino is disfavored. 

The situation is different if the gluino is heavy, rrig > 150 GeV. This gluino can be unstable due to weak R-parity 
violation |3(J] and have a lifetime Tg > 1 yr, i.e. long enough to be UHE carrier (see beginning of this section). 
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Then the calculated relic density at the time of decay is not in conflict with the cascade nucleosynthesis and all 
cosmologically produced g-hadrons decayed up to the present time. Moreover, the production of these gluinos by 
cosmic rays in the atmosphere is ineffective because of their large mass. 

Glueballino, or more generally <?-hadron, looses its energy while propagating from a source to the Earth. The 
dominant energy loss of the g-hadron is due to pion production in collisions with microwave photons. Pion production 
effectively starts at the same Lorentz-factor as in the case of the proton. This implies that the energy of the GZK 
cutoff is a factor m§ /m p higher than in case of the proton. The attenuation length also increases because the fraction 
of energy lost near the threshold is small, fi/nrig, where /i is a pion mass. Therefore, even for light <?-hadrons, 
irig > 2 GeV, the steepening of the spectrum is less pronounced than for protons. 

The spectrum of g-hadrons from the cosmic necklaces accounted for absorption in intergalactic space, is shown in 
Fig. 3. 

A very light UHE g-hadron interacts with the nucleons in the atmosphere similarly to UHE proton. The cross- 
section is reduced only due to the radius of g-hadron and is of order of ~ 1 rab [ff9| . In case of very heavy g-hadron 
the total cross-section can be of the same order of magnitude, but the cross-section with the large energy transfer, 
relevant for the detection in the atmosphere, is very small |ps)| . This is due to the fact that interaction of gluino in 
case of large energy transfer is characterized by large Q 2 and thus interaction is a deep inelastic QCD scattering. 

Thus, only UHE gluino from low-mass window 1 GeV < rrig < 4 GeV could be a candidate for observed UHE 
particles, but it is disfavored by the arguments given above. 



V. CONCLUSIONS 



Topological Defects naturally produce particles with extremely high energies, much in excess of what is presently ob- 
served. However, the fluxes from most known TD are too small. So far only necklaces (23| and monopole-antimonopole 
pairs |2^| can provide the observed flux of UHE CR. 

Another promising sources of UHE CR are relic superheavy particles [^4|,^5| . These particles should be clustering in 
the halo of our Galaxy |2f|, and thus UHE CR produced at their decays do not have the GZK cutoff. The signatures 
of this model are dominance of photons in the primary flux and Virgo cluster as a possible discrete source. 

Apart from protons, photons and neutrinos the UHE carriers can be neutralinos J28|, gluino p6|-p8[ and monopoles 
^9j,^7|. While neutralino is a natural candidate for the Lightest Supersymmetric Particle (LSP) in SUSY models, 
gluino can be considered as LSP only phenomenologically. LSP are naturally produced in the parton cascade at the 
decay of superheavy X-particles. In case of neutralino both fluxes and cross-sections for interaction is somewhat lower 
than for neutrino. In case of gluinos the fluxes are comparable with that of neutralinos, but cross-sections for the 
production of observed extensive air showers are large enough only for light gluinos. These are disfavored, especially 
if the charged gluebarino is lighter than the neutral one jscl . 
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FIG. 1. Predicted proton (p) and gamma-ray (7) fluxes from necklaces in comparison with experimental data. 
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FIG. 2: Predicted fluxes from decaying X-particles: nucleons (p,p,n,n) from 
the halo (curve extragalactic protons (curve Ip Xtr ), photons from the halo 

(curve Z^ al °), and neutrinos from the halo and the extragalactic space (curve 
Il ot ). 
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